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We present an investigation of magnetic ordering in the two-dimensional S — 1/2 quantum magnet 
Cu(Pz)2(C104)2 using specific heat and zero field muon-spin relaxation (/i + SR). The magnetic 
contribution to the specific heat is consistent with an exchange strength of 17.7(3) K. We find 
unambiguous evidence for a transition to a state of three-dimensional long range order below a 
critical temperature Tn = 4.21(1) K using ^i + SR even though there is no feature in the specific 
heat at that temperature. The absence of a specific heat anomaly at Tn is consistent with recent 
theoretical predictions. The ratio of Tn/J = 0.24 corresponds to a ratio of intralayer to interlayer 
exchange constants of \ J' / J\ = 6.8 x 10 -4 , indicative of excellent two-dimensional isolation. The 
scaled magnetic specific heat of [Cu(Pz)2(HF2)]BF4, a compound with an analogous structure, is 
essentially identical to that of Cu(Pz)2 (010.4)2 although both differ slightly from the predicted value 
for an ideal 2D 5 = 1/2 Heisenberg antiferromagnet. 



PACS numbers: 75.50.Xx, 75.50.Ee, 76.75.+i 

I. INTRODUCTION 

The properties of systems described by the S = 1/2 
two-dimensional square lattice quantum Heisenberg an- 
tiferromagnet (2DSLQHA) modelA remain some of the 
most pressing problems in condensed matter physics. 
The 2DSLQHA model is described by the Hamiltonian 

ff = ]Tjs,.s 3 , (i) 

where J is the nearest neighbour in-plane superex- 
change interaction. Although long range magnetic or- 
der (LRO) does not occur within the model above zero 
temperature^, layered materials that are well described 
by this model inevitably possess some degree of interlayer 
coupling (quantified by a coupling constant J') that leads 
to a crossover to a regime of three-dimensional (3D) LRO 
at a non-zero Neel temperature Tn. The identification of 
Tn and hence the ratio /cb^n / J allows the evaluation of 



the extent to which a material may be described by the 
2DSLQHA model. 

The measurement of Tn is often problematical in 
anisotropic spin systems due to the reduced ordered 
moment that typify these materials. In addition, the 
short-range order that exists in the two-dimensional (2D) 
planes of a layered material above a 3D transition reduces 
the effective number of degrees of freedom involved in the 
transition, diminishing the expected anomaly in specific 
heat£. Our recent study^ of the quasi one-dimensional 
(ID) 5 = 1/2 chain compound CuPz(NOa)2 (where Pz 
is pyrazine (C4H4N2)) has shown that implanted muons 
are uniquely sensitive to the presence of magnetic or- 
der in quasi-lD materials. In this paper we present 
an investigation of the magnetic properties of the two- 
dimensional analogue of CuPz(N03)2, namely copper 
pyrazine perchlorate (Cu(Pz)2(C104)2). This material 
has been the subject of several studies aimed at elu- 
cidating its properties^&i&^iii, although hitherto, no 
evidence of a magnetic transition has been found. We 
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have carried out detailed specific heat and muon-spin re- 
laxation measurements. The latter provide unambigu- 
ous evidence for a transition to a state of 3D LRO. 
As expected from recent theoretical studies of highly 
anisotropic spin systems, our specific heat measurements 
show no anomaly at the critical temperature and are 
shown to be in very good agreement with the recent pre- 
dictions of Monte Carlo simulations of the 2DSLQHA 3 . 

Copper pyrazine perchlorate (Cu(Pz) 2 (C10 4 )2) has 
long been held to be a good example of a 2DSLQHA 5 . 
The material is formed from layers consisting of rect- 
angular arrays of Cu 2+ ions bridged by pyrazine ligands, 
which provide the intralayer superexchange— . The 300 K 
structure of Cu(Pz) 2 (0104)2 has the Cu 2+ ions semico- 
ordinatcd with two disordered, crystallographically in- 
equivalent, CIO4 tetrahedra (one above and one below 
the layer) with C2/m symmetry^. The CIO4 ions or- 
der at low temperatures, where the crystal structure is 
reduced to C2/c, with two the sets of identical pyrazine 
pairs tilted by 62.8° and 69.1° with respect to the Cu- 
N coordination planed. This results in a herring-bone 
configuration, rather than a square lattice. Despite this 
distortion, magnetic susceptibility measurements fit very 
well to the 2DSLQHA mode£&±2 with J = 17.5 K and 
g = 2.11. 

A second copper pyrazine compound, 
[Cu(Pz)2(HF2)]BF4, has recently been reported^. 
X-ray measurements^ show it to have a similar layered 
structure to Cu(Pz) 2 (0104)2, but with adjacent layers 
bridged by HF^ ions. The intralayer exchange strength 
and 3D ordering temperature were found to be 5.7 K 
and 1.54 K, respectively. 



II. EXPERIMENTAL DETAILS 
A. Specific heat 

The heat capacity measurements on Cu(Pz)2 (0104)2 
have been carried out on four small crystals with total 
mass of 75.7 mg fixed on a Cu sample holder contain- 
ing a Constantan heater and a calibrated Cernox ther- 
mometer with Apiezon-N grease. We used a home-made 
calorimeter attached to a Janis He 3 refrigerator, with an 
automated system based on the semi-adiabatic method. 
The contribution of the sample holder as well as of the 
grease to the total heat capacity had been subtracted. 
Specific heat measurements for [Cu(Pz)2(HF2)]BF4 have 
been reported previouslyii. 

B. Muon-spin relaxation 

Zero field muon-spin relaxation (ZF /^ + SR) mea- 
surements have been made on a powder sample of 
Cu(Pz) 2 (C10 4 ) 2 using the MuSR instrument at the ISIS 
facility, Rutherford Appleton Laboratory, UK. In a /x + SR 
experiment— spin-polarized positive muons are stopped 



in a target sample, where the muon usually occupies an 
interstitial position in the crystal. The observed property 
in the experiment is the time evolution of the muon spin 
polarization, the behaviour of which depends on the lo- 
cal magnetic field at the muon site. Each muon decays, 
with an average lifetime of 2.2 /is, into two neutrinos 
and a positron, the latter particle being emitted prefer- 
entially along the instantaneous direction of the muon 
spin. Recording the time dependence of the positron 
emission directions therefore allows the determination of 
the spin-polarization of the ensemble of muons. In our 
experiments positrons are detected by detectors placed 
forward (F) and backward (B) of the initial muon po- 
larization direction. Histograms Ap(t) and -/Vb(£) record 
the number of positrons detected in the two detectors 
as a function of time following the muon implantation. 
The quantity of interest is the decay positron asymmetry 
function, defined as 

= N F (t) - a c , P N B (t) 
V ' N F (t) + a exp N B (t) V ' 

where a eX p is an experimental calibration constant. A(t) 
is proportional to the spin polarization of the muon en- 
semble. 

For these measurements, an Oxford Instruments Var- 
iox cryostat and a sorption pump cryostat were usecU 3 -, 
with a powder sample packed in a silver foil packet (Ag 
thickness 25 /j,m) and mounted on a silver backing plate. 
Silver is used since it possesses only a small nuclear mo- 
ment and so minimizes any background depolarizing sig- 
nal. 



III. RESULTS 
A. Specific heat 

The magnetic specific heat of the 2DSLQHA has been 
calculated using a variety of techniques. Early efforts 
concentrated on high temperature series expansions 14 
while more recent work has involved quantum Monte 
Carlo simulations 3 ^. These studies reveal that the spe- 
cific heat rises rapidly from zero at low temperature to 
a rounded maximum (of magnitude 3.8 J/K-mol) at a 
temperature T max of 0.65J before decreasing at higher 
temperatures. A rough approximation to the exchange 
strength is thus J = T max /0.65 = 1.54 T max . The pres- 
ence of interlayer interactions J' induces a transition to 
long-range order and causes the existence of a local max- 
imum in the specific heat at Tjy; this maximum is unde- 
tectably small for ratios \ J' / J\ < 0.01 but grows rapidly 
as the ratio approaches unity 3 -. 

The magnetic contribution of the ideal 2DSLQHA 
(J 7 = 0) can be represented as a function of tempera- 
ture as a ratio of polynomials in powers of T / J using the 
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FIG. 1: Specific heats of Cu(Pz) 2 (CI0 4 ) 2 (©) and 
[Cu(Pz)2(HF2)]BF4 (x) as functions of temperature. The 
solid lines correspond to the best-fit predictions for the spe- 
cific heats based on the parameters given in the text. The 
dashed lines represent the estimated lattice specific heats for 
Cu(Pz) 2 (C10 4 )2 (short dashes) and [Cu(Pz) 2 (HF 2 )]BF 4 (long 
dashes). The arrows mark the temperatures of the magnetic 
ordering transitions for the two compounds as determined by 
muon spin relaxation studies. Note the absence of anomalies 
in the specific heats at those temperatures. 



for this compound 11 -. The data sets for both com- 
pounds have been analyzed by assuming the specific 
heats have both lattice and magnetic contributions. The 
lattice contribution in the low temperature region can 
be approximated as Ci a tticc=aT 3 + PT 5 + jT 7 , where 
a, 0, 7 are constants to be determined. By model- 
ing the Cu(Pz)2 (0104)2 data over the entire range as a 
sum of the magnetic and lattice contributions, we ob- 
tain the best fit for the parameters J = 17.7(3) K, 
a = 0.0039(2) J/K 4 -mol, /? = -2.75xl0~ 6 J/K 6 -mol, 
and 7 = 6.38xl0~ 10 J/K 8 -mol, while the correspond- 
ing parameters for [Cu(Pz)2(HF2)]BF4 are J = 5.6(1) K, 
a = 0.0114(5) J/K 4 -mol, (3 = -6(1) xl0~ 5 J/K 6 -mol, 
and 7 = 1.2(2) xl0~ 7 J/K 8 -mol. The resultant fits to 
the two compounds are shown as the solid lines while the 
estimated lattice contributions for Cu(Pz)2 (0104)2 and 
[Cu(Pz)2(HF2)]BF4 appear as the lines of short dashes 
and long dashes, respectively. The value of the exchange 
strength J is in excellent agreement with the value of 
17.5 K obtained from analysis of the magnetic suscep- 
tibility data£ i 8 ' 9 ' 10 and from the dispersion relation as 
obtained by inelastic neutron scattering^. Likewise the 
value of 5.6 K obtained for [Cu(Pz)2(HF2)]BF4 agrees 
well with the value of 5.7 K obtained from susceptiility 
studies^. 



following equation 
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where R is the gas constant and C n , D n are coefficients 
obtained by fitting data from Monte Carlo Simulation 3 . 
The coefficients appear in Table HI This form is useful 
for comparing experimental data to the theoretical pre- 
dictions and is valid for T/J > 0.1. 
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( h 


D n 


1 


-0.0036 


1.86131 


2 


0.26197 


-11.51455 


3 


-1.45326 


28.45991 


4 


3.06373 


-32.81602 


5 


-0.38349 


18.6304 



TABLE I: Fitting coefficients for pure 2DSLQHAF specific 
heat polynomials from Ref. y. 



The molar specific heat data for Cu(Pz)2 (0104)2, rep- 
resented as circles, are shown as a function of temper- 
ature in Fig. [TJ There is no low-temperature maxi- 
mum in the data but an inflection point occurs near 
8 K, indicative of a non-lattice contribution. In contrast, 
the specific heat of the analogous 2D copper pyrazine 
compound, [Cu(Pz) 2 (HF 2 )]BF4, also shown in Fig. [TJ 
shows a definite low-temperature anomaly, correspond- 
ing to the lower exchange strength of 5.7 K reported 



The magnetic specific heats, obtained by subtracting 
estimated lattice contributions from the total specific 
heats of Fig.[TJ are shown in Fig. 2(a)| along with the es- 
timated lattice contributions as the dashed lines, plus the 
best fit to the theoretical prediction for the 2DSLQHAF. 
While the overall agreement between data and theory 
is quite good, there are systematic disagreements that 
are similar for each compound. The theoretical predic- 
tion initially rises more slowly at low temperature, sur- 
passes the experimental data while reaching its maxi- 
mum value at a lower temperature, then decreases less 
rapidly than the data at higher temperatures. Neverthe- 
less, the entropy changes for both data setsil and the 
theoretical curve are all within several percent of i? ln2. 
It is important to note that neither data set shows an 
anomaly at the temperature of the 3D ordering temper- 
atures (4.21 K for Cu(Pz) 2 (C10 4 )2 [see below], 1.54 K 
for [Cu(Pz) 2 (HF 2 )]BF 4 [see Ref. Gt[) as determined by 
At+SR. 

The comparison of experiment and theory is more 
clearly represented in Fig. 2(b) in which the reduced spe- 
cific heat (C mag / R) is plotted as a function of the reduced 
temperature (T/J) for the theory and both experimental 
data sets. Here it is noticed that the two experimental 
data sets are essentially supcrimposible up to T « J. 
At higher temperatures, the data for Cu(Pz) 2 (0104)2 
falls below that of [Cu(Pz) 2 (HF 2 )]BF 4 and of the theory. 
However, by this temperature, the lattice contribution 
of Cu(Pz)2 (0104)2 is more than 90% of the measured 
value, with commensurate uncertainties in the magnitude 
of C mae: . 
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FIG. 2: (a) Magnetic specific heats of Cu(Pz) 2 (CI0 4 )2 (0) 
and [Cu(Pz)2(HF2)]BF4 (x) as functions of temperature with 
the solid lines corresponding to the best-fit predictions for the 
magnetic specific heats based on the parameters given in the 
text. The dashed lines (short dashes for Cu(Pz)2 (0104)2 and 
long dashes for [Cu(Pz)2(HF2)]BF4) mark the estimated lat- 
tice contributions that were subtracted from the data sets in 
Fig.[T]to produce the magnetic specific heats, (b) The reduced 
magnetic specific heats (C m&g /R) of Cu(Pz) 2 (C104)2 (©) and 
[Cu(Pz)2(HF2)]BF4 (x) as functions of the reduced tempera- 
ture T/J with the solid line corresponding to the theoretical 
prediction 3 for the magnetic specific heat of the 2DSLQHA. 



B. n+SK 

ZF /i + SR spectra measured on Cu(Pz)2 (0104)2 at two 
temperatures are shown in Fig. [3l For temperatures 
T > 4.27 K the measured spectra are found to contain 
two contributions. The first is a fast relaxing compo- 
nent Ai which dominates the signal at early times and is 
well described by an exponential function exp(— Xt). The 
second is a larger, slowly relaxing component A2, which 
dominates at intermediate times and fits to the Kubo 
Toyabe (KT) function^ /kt(A,4) where A is the second 
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FIG. 3: (a) ZF /i+SR spectra in Cu(Pz) 2 (C10 4 )2 measured 
at T = 0.35 K and 5.02 K. At temperatures T > 4.27 K, Eq. 
(UJ describes the data while the oscillations measured below 
this temperature are fitted to Eq. ([5j|. 

moment of the static, local magnetic field distribution 
defined by A = 7 M \/((B - (_B)) 2 ), where B is the magni- 
tude of the local magnetic field and 7 M (= 27r x 135.5 MHz 
T -1 ) is the muon gyromagnetic ratio. The KT function is 
characteristic of spin relaxation due to a random, quasi- 
static distribution of local magnetic fields at diamagnetic 
muon sites. We do not observe the recovery in asymmetry 
at late times that is expected for the static KT function. 
The lack of this recovery is probably due to slow dynam- 
ics in the random field distribution and is crudely mod- 
elled here with an exponential terrn^ exp(— K^t). The 
data were found to be best fitted with the resulting func- 
tion 

A(t) = Aiexp(-Aii) + A 2 / KT (A, i) exp(-A 2 t) 

+A bg cxp(-Abg), (4) 

where Abg exp(— Xbgt) represents a constant background 
signal from those muons that stop in the sample holder 
and cryostat tail. The fitting parameters show very little 
variation in the range 4.27 < T < 14 K. The compo- 
nent with amplitude A\ (found to have relaxation rate 
Ai ~ 10 MHz) very likely arises due the existence of a 
paramagnetic muon state, fn the component with ampli- 
tude A2, the magnitude of A(~ 0.2 MHz) suggests that 
the random magnetic field distribution giving rise to the 
KT function is due to nuclear magnetic moments, im- 
plying that the field due to electronic moments at these 
muon sites is motionally narrowed out of the spectrum 
due to very rapid fluctuations. 

In spectra measured at temperatures T < 4.20 K oscil- 
lations in the asymmetry spectra are observed at several 
frequencies (Fig. [3]) . These oscillations are characteristic 
of a quasistatic local magnetic field at the muon stopping 
site, which causes a coherent precession of the spins of 
those muons with a component of their spin polarization 
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i Aj (%) <fo 

1 3.127 -39.5 

2 2.508 -21.88 

3 0.174 96.6 

4 2.633 



TABLE II: Fitting parameters for equation <(Sj applied to data 
measured for T < 4.20 K. 



perpendicular to this local field (expected to be 2/3 of 
the total polarization). The frequency of the oscillations 
is given by Vi — 7^i?.;/27r, where Bi is the magnitude 
of the local magnetic field at the ?th muon site. Any 
distribution in magnitude of these local fields will result 
in a relaxation of the oscillating signal, described by re- 
laxation ratesi^ A^. The presence of oscillations at low 
temperatures in Cu(Pz) 2 (0104)2 suggests very strongly 
that this material is magnetically ordered below 4.20 K. 

Three separate frequencies were identified in the low 
temperature spectra, corresponding to three magneti- 
cally inequivalent muon sites in the material. The preces- 
sion frequencies, which are proportional to the internal 
magnetic field as experienced by the muon, are propor- 
tional to the magnetic order parameter for these system. 
In order to extract the T-dependence of the frequencies, 
the low temperature data were fitted to the functional 
form 

A(t) = Ai exp(— Ait) cos(27Wii + 4>i) 
+A 2 cxp(-A 2 i) cos(2wv 2 t + &) 
+A3 cxp(-A 3 i) cos(2iris 3 t + cf> 3 ) 
+A 4 cxp(-A 4 i) + A hg exp(-Abgi), (5) 

We note that phase offsets (j> were required to fit the data 
(see Table HH as observed in our previous ^ + SR studies 
of Cu 2+ system o 4 ' 20 i 21 . The amplitudes of the oscillat- 
ing components were found to be constant in the ordered 
state and were fixed at the values shown in Table [UJ 
These values suggest that the three muon sites have oc- 
cupation probabilities in the ratio 0.60 : 0.33 : 0.07. The 
term A4 exp(— X^t) accounts for the contribution from 
those muons with a spin component parallel to the lo- 
cal magnetic field expected to be half of the oscillat- 
ing amplitude (see above). From Table HTl we see that 
A±/{Ai + A2 + A3) w 0.45 (i.e. the ratio of amplitudes 
resulting from local magnetic field components parallel 
to the initial muon-spin direction to those perpendicu- 
lar), is close to the expected value of 1/2, suggesting that 
the material is ordered throughout its bulk. Further ev- 
idence for the presence of a magnetic phase transition is 
provided by the observation that the relaxation rates Ai 
tend to increase as Tn is approached from below, due to 
the onset of critical fluctuations. 

The three frequencies were found to be in the propor- 
tions v\ : vi : vj, = 1 : 0.56 : 0.12 and were fixed in 
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FIG. 4: Temperature evolution of the three magnetic os- 
cillation frequencies in Cu(Pz)2(C104)2, extracted from fits 
to Eq.©. The solid lines are fits to the functional form 
Vi(T) = !/i(0)(l - (T/T N ) Q )' 3 (see main text). 

this ratio in the fitting procedure. The magnitude of 
these frequencies were fitted as a function of tempera- 
ture with parameters in Eq.((5]) fixed at the values given 
in Table |TTJ and relaxation rated allowed to vary. The re- 
sulting temperature evolution of the precession frequen- 
cies is shown in Fig. [4] From fits of the data to the 
form Vi{T) = 1^(0) (1 - (T/T N ) Q )' 3 , we estimate T N = 
4.21(1) K, a = 1.8(3), = 0.29(2), z^(0) = 2.38(3) MHz, 
iaj(O) = 1.33(2) MHz and 1/3 (0) = 0.29(2) MHz. 

IV. DISCUSSION 

A method for estimating the interlayer coupling con- 
stant J' in 3D arrays of 2DSLQHA has recently been de- 
veloped, based on a modified random phase approxima- 
tion, modelled with classical and quantum Monte Carlo 
simulations 22 . This approach leads to an empirical for- 
mula relating J' and Tn, 

|jyj|=exp(&-^), (6) 

where p is the spin stiffness, given 2 ^ by p s = 0.183 J 
and b = 2 A3 for S = 1/2. Using J = 17.7 K, we ob- 
tain p s = 3.26 K and \ J'/J\ = 6.8 x 10~ 4 . Table [En] 
compares these parameters to those for other 2D layered 
compounds. We see that Cu(Pz)2(C104)2 compares very 
favourably as a highly anisotropic system, with the lay- 
ers around an order of magnitude better isolated than 
those in Cao.ssSro.isCuO^, although not as successful a 
realisation of a 2DSLQHA as Sr 2 Cu0 2 Cl2, where Eq.© 
yields \ J' / J\ ~ 10~ 5 .The J' /J ratios were obtained as- 
suming perfectly Heisenberg exchange interactions. The 
presence of any exchange anisotropy will also raise the 
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Parameters for 
Cu(Pz) 2 (C10 4 )z 



TABLE III: 
SraCuOaCl^ 

Cao.ssSro.isCuO^, '(5MAP) 2 CuBr4^ 
estimated from Eq.©. 



the layered compounds 
9 rCu(Pz) 2 (HF 2 )]BF*ii, 
with values of \J' / J\ 



\J\/ka (K) T N (K) \k B T N /J\ \J'/J\ 



Sr 2 Cu0 2 Cl 2 1451 
Cu(Pz) 2 (C10 4 ) 2 17.7 
[Cu(Pz) 2 (HF 2 )]BF 4 5.7 
Cao.85Sro.i5Cu0 2 1535 
(5MAP) 2 CuBr 4 6.5 



256.5 0.18 

4.21 0.24 

1.54 0.27 

537 0.35 

3.8 0.58 



2.5 x 10 _t> 
6.8 x 10" 4 
2.3 x 10" 3 

1.6 x 10~ 2 
2.1 x 10" 1 



ratio Tn / J so by ignoring the ansiotropy fields while us- 
ing Eq.([6]), we are setting an upper bound for the J' / J 
ratios. 

Cu(Pz) 2 (C10 4 )2 and [Cu(Pz) 2 (HF 2 )]BF 4 contain 
structurally similar copper pyrazine layers separated ei- 
ther by the bulky perchlorate groups or by the HF 2 an- 
ions. Their intralayer exchange strengths differ by a fac- 
tor of three (17.7 K, 5.6 K) but their low ratios of T/v/J 
show them to be well isolated (J' /J ~ 10~ 3 , see Ta- 
ble As seen in Fig. |2(b)[ their magnetic specific 
heats are essentially identical up to relative tempera- 
tures of T » J. Nevertheless, there is noticable differ- 
ence between their behavior and that predicted for iso- 
lated 2DSLQHA model systems^. This difference may 
arise from the small terms in the Hamiltonian relevant 
for the compounds but ignored in the Monte Carlo sim- 
ulation. Studies of Cu(Pz) 2 (C104) 2 in the magnetically 
ordered state show the presence of about a 0.3 T axial 
field perpendicular to the magnetic layers^. While small 
compared to the 60 T exchange field, this axial field will 
break the rotational symmetry of the Heisenberg Hamil- 
tonian and will modify the low-temperature behavior. It 
will be necessary to examine the specific heat of the quasi- 
2DSLQHA in the presence of both weak interlayer inter- 
actions as well as weak anisotropics to understand these 
effects qualitatively. 



V. CONCLUSION 

We have examined the 2D Heisenberg antiferromag- 
net Cu(Pz) 2 (C104) 2 using specific heat and /x + SR tech- 
niques. The value of the exchange strength obtained (J — 
17.7(3) K) from the magnetic specific heat is in excellent 
agreement with values obtained from susceptibilitjA&SiiS 



and neutron scattering^ experiments. The muon spin 
relaxation studies demonstrated the existence of 3D 
long range magnetic order throughout the bulk of the 
sample in Cu(Pz) 2 (C104) 2 at temperatures below Tn 
= 4.21(1) K. The value of T N /J (0.24) corresponds 
to a very low ratio of interlayer to intralayer ex- 
change strengths (| J' / J\=6.8x 10 -4 ). The combination 
of good isolation of the magnetic layers plus relatively 
small exchange strength makes Cu(Pz) 2 (C104) 2 , as well 
as [Cu(Pz) 2 (HF 2 )]BF4, good candidates for studies of 
the field dependence of the energy spectrum of the 
2DSLQHA2&22. 

The magnetic ordering transition is readily apparent 
in the /i + SR data (Fig. [3]) even though no anomaly is 
apparent in the specific heat results at that temperature 
or at T N of [Cu(Pz) 2 (HF 2 )]BF 4 (Fig.[2(a)j). The absence 



of an ordering peak in the specific heat for well isolated 
low-dimensional magnets was previously anticipated 3 - in 
Monte Carlo studies. These demonstrate that most of 
the available spin entropy (i?ln2 per mole) is associ- 
ated with the development of in-plane correlations which 
give rise to the broad peak in C mag centred around 
T = 0.65J. This leaves little remaining entropy to be 
associated the 3D ordering at T = Tn- The critical 
temperature of Cu(Pz) 2 (C104) 2 was also not detected 
by previous inelastic neutron scattering experiments on 
polycrystalline samples; only after adequate single crys- 
tals became available could the transition be observed 
with this technique 16 . These results emphasize the value 
of muons in detecting magnetic order in polycrystalline 
materials. 
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